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I.-- INTRODUCTION 


The following report is a summary of the design configuration and 
performance characteristics of a Saturator designed and constructed for the 
NASA. The hardware was developed to provide ground-based simulation for 
some of the experiments for ACPL-1 first flights of Spacelab. The design 
shared concepts and technology developed somewhat earlier for the Continu- 
ous Flow Diffusion Chamber (CFD)[5], Some of the difficulties encountered 
with the apparatus during this project are discussed, and recommendati ons 
concerning testing of this type of instrument are presented in the text as 
well as in Appendix I. A theoretical study, undertaken before construction 
of the instrument, is presented in Appendix II. (Part of Appendix II 
appeared in earlier progress reports of this Contract). 

Two Saturators were constructed; one instrument was integrated into 
a common control system with the CFD, and the second unit was delivered to 
the NASA with a comparatively simple manual control system. Some minor 
improvements were incorporated into the second Saturator based upon experi- 
ence with the first device. 

The new Saturators provide a means of accurately [3] fixing the water 
vapor mixing ratio of an aerosol sample. Dew point temperatures from almost 
freezing to ambient room temperatures can be attained with high precision. 

The instruments can accommodate aerosol flow rates approaching 1000 cc/s. 
Provisions were made to inject aerosols upstream of these Saturators* although 
downstream injection can be accomplished as well. 

A device of this type will be used in the ACPL-1 to condition various 
aerosols delivered concurrently to a CFD, Expansion Chamber, and Static Diffu- 
sion Chamber used in zero gravity cloud-forming experiments. The Saturator 
was designed to meet the requirements projected for the flight instrument as 
specified in SVS-9800 [3], 
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n. DESCRIPTION OF THE INSTRUMENT 

The Saturator incorporates a flow passage consisting of two closely 
spaced parallel plates maintained at a stable and well-known temperature. 

The plates are maintained wet; water is. supplied to the plates through a 
capillary duct connecting the plates to a water storage sump. The length 
of the flow passage, is such that approximately ten time constants' are: avail — 
aDie tor cnrtusion of water vapor into the moving airstream at flow rates 
less than about. 1000 cc/s. 

Temperature control of the plates is accomplished by circulating a 
coolant through channels machined in the thermal plates. This coolant is 
separated from the chamber space by appropriate seals. The thermal plate is 
constructed of aluminum plate; the coolant consists of a 4 to 1 mixture 
of distilled water and Prestone antifreeze. Dye present in the antifreeze 
solution provides a reliable means of determining if coolant leaks exist 

ID tfl6 arr, ^ ov> ° n3 ra 

ATst'agnant. air region approximately 2^5 cm wide is provided, around 
both longitudinal edges of the chamber space in order to isolate the aerosol 
from thermal edge effects. Ample insulation is provided over the thermal 
plates to minimize thermal losses to the environment. Chamber dimensions are 
shown in Figure 1. A plastic body fitted with 0.13 mm diameter stainless 
steel wires stretched; transversely across the channel of flow- is mounted just 
•downstream of the wetted surfaces of the thermal plates. These wires serve 
as electrical heater elements and are provided for the purpose of raising 
the temperature of the sample aerosol to prevent condensation downstream 
of the Saturator. The heater body is constructed of thermoplastic (ABS) with 
a low thermal conductivity in order to restrict heat conduction from the 
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CHAMBER DIMENSIONS 


Plate Width 38 cm 
Chamber Space Width 33 cm 
Plate Spacing 0.15 cm 
Wetted ‘Length 27 cm 


Post Heater Body Length 7.0 cm 


ro 

i 

ro 


LEGEND : 

1. Air Inlet Manifold 

2. Thermal Plate 



3. Stainless Steel Wicking Surface 

4. Copper Coated ABS Post' Heater Body 

5. Post Heater Elements 


6. Exhaust Manifold- 


Figure 1. Cross-Sectional Schematic of Saturator 



post heater body to the thermal plates. Radiation heat transfer from 
the wires to the thermal plates is not significant because the transverse 
wires operate at just a few degrees above the temperature of the Saturator, 
and because of the small surface area of the wires. The ABS post heater 
body is electroplated with a thin layer of copper in order to prevent 
electrostatic charging of these surfaces. 

Large metal manifolds are fitted to the inlet and outlet of the" 
Saturator to provide a plenum effect for the air entering and leaving the 
Saturator. Aerosol is introduced using turbulent mixing at the entrance 
of the upstream manifold. A pressure port is fitted in the exit manifold. 

Temperature measurements are taken on the back side of one thermal 
plate at downstream locations near the coolant inlet manifold. Holes 
are milled in the heavy aluminum thermal plate to accept closely fitting 
brass holders coated with thermal grease. Glass bead thermistors were 
potted into the holders using thermally conductive epoxy. During calibra- 
tion testing of the. instrument, a larger hole was machined into the thermal 
plate to accept a quartz thermometer probe, as discussed in Appendix I. 


2*?3 



o?X' 


C 


pffcj&t 


***** 


not 




111 • TEMPERATURE CONTROL 

As mentioned earlier, the first Saturator shares the control sys- 
tem developed concurrently for the CFD. The system is shown in Figure 2. 
Control hardware incorporated an INTEL 8010 microprocessor. 

A flow schematic for the first Saturator is shown in Figure 3. Cool- 
ant is circulated through the thermal plates and a thermoelectric powered 
heat exchanger. .Circulation is provided by a small centrifugal- pump at a 
rate of about 30 cm /s to each plate. Just upstream of the circulation pump 
a small resistance Servo heater is fitted. The coolant is also circulated 
over a small glass 'bead thermistor just downstream of the pump; this method 
of temperature measurement provides an accurate estimate of the bulk coolant 
temperature because of the thorough mixing provided by the pump. The sensi- 
tivity and short time constant available when using the thermistor provides 
an accurate measurement of the rate of change of the bulk coolant tempera- 
ture. This information is used to provide feedback for the immersion heater 
circuitry. Temperature stability better than ±0.01°C has been achieved-' 
using- this system. More detailed discussions of temperature measurement 
performance and temperature calibration with these new systems are provided 
in Reference 5.. 

The thermoelectric powered heat exchangers originally failed to per- 
form satisfactorily because of inadequate clamping between the heat exchanger 
and TEM's; this problem was solved by redesigning the clamping method. Addi- 
tional cooling performance was obtained by reducing the hot side tempera- 
ture of the TEM's. These problems and solutions are discussed in more detail 
in Reference 5. 

The second Saturator is provided with a Haake Model FK laboratory 
circulator instead of the system described above. 
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Figure 2a. Photo of Saturator 
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Figure 2b. Front View of CFD and 
Saturator control rack 


Air flow control panel 
Royco 225 mainframe 

Computer and power supply 

28 VDC supply 

TEM waste heat radiators 


Teletype 


3-2 


preceding page blank not filmed 

IV. WETTABLE METAL SURFACES 

The wettable metal surfaces used in the Saturator consist of a 
random, diffusionally. bonded mesh .of micronic-sized stainless steel fibers. 

The media is 0.32 mm thick and is purchased as Type 304 stainless, steel. 

Special Media No. 86 from Fluid Dynamics of DeLand, Florida. The media 
is. subjected to a special- Fluid .Dynamics process to reduce carbon contact 
to 0.04%. 

The media is diffusionally bonded to a stainless steel sheet .in a 
hydrogen-fired oven at a temperature of 980 to 1090°C. During the bonding 
process, the stainless media is in contact with the substrate and is under 
some pressure to assure uniform adhesion. In order to reduce corrosion in 
the finished- wi eking surface and substrate, the part undergoes a second 
heat treatment after bonding;: in a hydrogen oven fixtured to allow rapid 
cooling, the temperature of the part is raised to 1090°C and then rapidly 
quenched in N2»such that the time of exposure in the temperaature range of 
425 to 81 5° C is below about. 5 minutes. 

Part distortion was experienced, using an 0.8' mm thick substrate. 
Substrates 1.6 mm thick were acceptably flat, after processing. ' Care must be 
exercised during the diffusion bondina orocess not to stretch the substrate. 

These surfaces perform extremely well in the Saturator and CFD instru- 
ments. Uood wettability is achieved, with- the- supply sumps located at- the 
bottom of the instrument" in 1 g. During the Saturator evaluation testing', 
all tests were fun with the sump in the “down" position. Furthermore, dry 
air was .run through the instrument for a period of about one hour before 
taking data to make sure, the rate of flow of distilled water to the wicking 
surfaces was adequate. 
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A problem has been experienced with the sump design just after 
filling that is not fully explained; just after filling (perhaps over- 
filling) the sump, some liquid water can sometimes be drained from the 
Saturator (a few cm ). After the initial draining, the problem does not 
recur. 

At an air flow of 500 cc/sec, the Saturator can be run over four 
hours before it is necessary to refill the sump. The sump is packed with 
spun glass and has dimensions of 7.1 cm diameter x 30 cm length. 

The wicking surfaces are cleaned using the glow discharge cleaning 
process. Based upon experiences with the CFD and Saturator, it appears 
that several months of use can be obtained from these instruments before 
cleaning is necessary if adequate filtering is present. Our filters con- 
sist of activated charcoal and were designed to remove any long chain hydro- 
carbons. 

Handling of the wicking surfaces is kept to a minimum during cleaning 
and assembly. All leak testing is performed before cleaning, and assembly 
of the unit is completed as far as possible before the start of the cleaning 
operation. Using these concepts, satisfactorily clean surfaces can be 
obtained [5], 
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V.. CORROSrON PROBLEMS 

Initially the post heater body was plated with chrome over a base of 
nickel. Severe corrosion took place on this surface during use in the Satura- 
tor. This surface was stripped off on the first unit. The post heater body 
-on both instruments is now coated only with a thin layer of copper. This 
seems to have eliminated the corrosion problem encountered with the first 
unit. 

uorrosion pitting of the aluminum thermal plate was eliminated 
by using a 4/1 mixture of distilled water and Prestone antifreeze. Alumi- 
num and copper samples soaking in this mixture for over a period of a year 
) 

have shown no sign of corrosion. No pitting has been seen on the thermal 
plates. ■ Distilled water alone will cause severe pitting on aluminum in 
just a short period of time if any heavy metal ions are in solution. 

Some pitting has been experienced where aluminum is in close con- 
tact with stainless steel -and is exposed to the distilled water supplied to 
the wickincj surfaces. - . RirkT fjlIiAH 5 

PRE cE D -naP^^ N0TF 
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APPENDIX I 


SATURATOR CALIBRATION BY 
A GRAVIMETRIC METHOD 
Dr. F. Rogers 

I. PHILOSOPHY 

A measured mass of dry air was passed through the saturator at; a flow 
rate of approximately 500 cm 3 sec~^ , while the output -of the saturator .passed 
through sealed dessicant canisters. The canisters were weighed before and 
after each experiment to give the weight of trapped water. Two simple 
checks were routinely utilized to determine whether or not the dessicant 
traps were completely efficient. The measured mixing ratio was - then compared' 
with theory 

II. METHOD 

A. Weighing of Dry Air 

It was origin ally intended to complete all tests with dry air supplied 

i - 

by the Matheson Company and certified to have a dew point temperature no higher 
than “60°C. Tests 1 through. 4 were conducted wtth this air, but it was later 
necessary to use locally-available, dried, compressed air when the decision 
was made to continue with tests 5 through 8. In the former case^ the amount of 
water vapor present in the gas presents an estimable but negligible error in 
the measurement of cue- mass or ary air passed tnrougn the saturator. In the 
latter case, it was necessary to measure the water vapor content of the 
locally-obtained air. 

For all experiments, the air was supplied from pressure bottles weighing 
about 13.7 Kg. An additional 5.4 Kg was added by the pressure regulator 
and connections necessary to interface between the bottle and the saturator. 
Fortunately, it was possible to obtain an excellent high-capacity balance. 



the Ohaus Model 1119D Heavy Duty Solution Balance, of total capacity 20 kg 
and- sensitivity 1 g. The sensitivity was verified by loading the balance 
with a typical air bottle, then adding a 1 g weight after the balance had 
been zeroed; a detectable deflection was always found to result. In actual 
use, the balance, loaded with air bottle plus regulator, was zeroed before 
each experiment. At the end of each run, masses calibrated on Mettler balan- 
ces were added to the pan to re-establish the zero - the poises supplied with 
the balance were not touched - and the mass of air expended was given by the 
calibrated masses. Therefore, the mass of air expended, usually around 900 g, 
was considered to be measured with an error equal to the sensitivity of the 
balance, or about + 1 g (the repeatability of the balance seemed to be better 
then +1 g). Care was taken to ensure that the hose connecting the air bottle 
to the Saturator did not cause spurious weight changes because of movement 
of the hose. Photo #1 shows the Ohaus balance zeroed with an air bottle in 
position on the pan, while mounted in a wooden cradle. 

The weight of water vapor contained in the gas must be subtracted 
from the weight loss of the cylinder to give the weight of the dry air. In 
the case of runs 1 through 4, using Matheson "dry air" of dewpoint <-60°C, 
the maximum error caused in the. measured mixing ratio would be less than 0.1%, 
since the saturation vapor density at -60°C is about .025 gms/m 3 . In the 
case of runs 5 through 8, the amount of water vapor present in the locally- 
supplied "dry air" was measured, using a bag .sample and one of the canisters 
described in Section G.l. The mixing- ratio, -to an -accuracy of only +80%, 
was found to be 0.5 g/kg. However, the resulting error in the dry air mass 
measurement was only about 0.04%, as typically about 900 g dry air were 
expended. In runs 5, 6 and 7, canisters of dessicant were inserted between 
the air source and the Saturator to intercept this moisture and allow its 
weight to be subtracted from the weight of the air expended, but this attempt 
was abandoned in run 8 due to continuing problems with leaks. 



B. Saturator Support Systems 
1. Temperature Control 

For all the runs (1 through 8), the Haake thermos tatted circulating 
bath, to be supplied with the unit, was used for plate temperature control, 
and was left at the same setting, arbitrarily about 19.6°C. In all cases, 
fluctuations greater than about +0.01 °C were not observed. 

^2. Temperature Measurement 

A special well was drilled in one of the saturator plates to receive 
the probe from a Hewlett-Packard Model 2801 A Quartz Thermometer, with abso- 
lute accuracy +0.05°C and resolution 0.001°C (see Photo #2). In the vici- 
nity of 20°C, and at about 900 mb total pressure, 0.05°C corresponds to an 
error of about 0.3% in mixing ratio. 

3, Pressure Measurement 

The absolute pressure inside the saturator was measured with a Setra 
Systems, Inc., Model 204 pressure transducer. This device operates over the 
range zero to .1700 mb, with the major error contributions, according to 
factory specs, being non-linearity (0.1%) and hysteresis (0.05%). The ori- 
ginal factory calibration graph was obtained, and showed that, for this par- 
ticular instrument, hvsteres.is was limited to about 0.02% in the operating 

region of interest. A far worse error was observed at this laboratory; the 
zero of the device, checked during and after the calibration runs,, seemed to 
drift by about 2 mb, or about 0.2% of the typical value measured in the 
saturator. Therefore, the pressure measurement error, estimated by the 
"root sum-of-squares" of the non-linearity, hysteresis, and zero-drift 
errors, was about 0.2%. This error contributes to the error of calculated 
mixing ratio, as discussed in Section III. A. 



4. Water Supply to Wi eking Surfaces 

The saturator sump was filled with distilled and deionized water. 
After run #4, it was found that the sump was nearly dry, after an estimated 
five hours of saturator operation at 19.6°C on the dry air discussed. 

5. Saturator Warm-up and Flush 

Starting with run #3, the practice was adopted of allowing ordinary, 
compressed, "breathing quality" aiir (relative humidity about 15%) to cir- 
culate through the saturator for at least one hour before the actual cali- 
bration run began. In runs 4 through 8, no more than two' minutes elapsed 
between shutting off this "flushing" flow and beginning the calibration flow. 
The temperature control „and measurement devices were always turned on at 
least two hours before calibration began. 

C. Measurement of Water Trapped Downstream of Saturator 
1. Dessicant Canisters 

Three 1000 ml "Nalgene" bottles, arranged in parallel, received 
the saturator output. These were filled with indicating "Drierite", and 
efficient dessicant which, in equilibrium, reduces the dewpoint of air to 
approximately -80°C. These canisters (see Photo #3 and #4) were fitted 
with inlet and outlet tygon tubing connections, manifolds to distribute 
incoming air into the Drierite packs, and were leak- tested several times 
at an overpressure of about 2 psi. The flow resistance of these canisters 

and of a flow meter also placed downstream of the saturator (on the canister 
exhaust) caused the saturator pressure to be elevated about 1 psi above 
ambient. The saturator was previously leak tested at about 5 psi. These 
canisters were separately weighed before and after each run on an ordinary 
2-pan laboratory balance (Ohaus Model 1600) sensitive to about 0.03 g. „ 

Masses calibrated on a Mettler H51AR balance were used for the measure- 
ment, but it was not possible to use the Mettler balance for the calibra- 
tion itself. The calibration procedure would improve considerably if all 



three canisters could be weighed at the same time, say on a top-loading 
balance of better sensitivty. 

Let 8-j be the weight of canister #1 before a calibration run, and A-| 
be its weight afterward. Then, the weight gain of the canister is given 
by G 1 

B 1 + 0.03 g 

A-|’± 0.03 g 

G-j ± >/o.03^ + 0.03^ g, 

or, the propagated error in G-j is about +0.04 g. The total water trapped by 
all three canisters is given by 

| Gi + /l (0.04) 2 , 

or the total propagated error is about + 0.07 g. Since about 14 g water were 
typically measured, the total error contribution to the measured mixing ratio- 
is- about to. 5%. 

2. Checks of Canister Efficiency 

Loose Drierite in a container can be, subject to "channeling"; that is, 
air can sometimes find preferred channels of relatively low flow resistance 
rather than evenly permeating the- entire dessicant pack.. Indicating Drierite 
changes from blue to red in color when water absorption, has taken Diace. and 
the first' indication that channeling had occurred, or that the amount of dessi- 
cant in a given canister was too small, would be a red zone penetrating to the 
outlet end of the canister. This visual check of canister efficiency always 
indicated that, in the first approximation, air coming from the saturator at 
about 500 cm sec had more than adequate residence time in the dessicant 
canisters for equilibrium with the Drierite. As a second check, however, a 



bag sample was taken downstream of the canisters "(and of a flow meter used 
for rough monitoring of the system air flow rate). The bag,, of about 200 L, 
capacity, could usually be "on-line" for about 8 minutes out of the total 
run length of about 30 minutes; it was not attached until each calibration 
run had been in progress about 5' minutes. Following the run, about two-thirds 
to. three-quarters of the bag contents could effectively be delivered to one 
of the Drierite. -canisters, .to give a measure of. how much, if any,, water had 
passed the three primary traps. Again, masses of water less than about .04 g 
could escape detection by this procedure (see Section C.l). Furthermore, since 
the bag procedure represented only about (2/3 x 8 =) 5 minutes out of the total 
run time of around 30 minutes, it was necessary to multiply the moisture con- 
tent of the bag by a factor of about 6. Therefore, if .04 g were present in 
the bag, the total moisture escaping the traps would be estimated to be 
around 0.2 g, or some 1.5% of the- total water collected in any one run. There- 
fore, this procedure represents the largest potential error source of the 
entire calibration effort - but the adjective "potential" is used because it 
seems unlikely, in the light of the first check mentioned above, that 0.2 g 
escaped detection every time. In future efforts, however, it is certainly 
desirable- to improve this check-orocedure. for example, ‘by utilizing a bag of 
Targe enough capacity to collect the exhaust air during the entire time of 
an experiment - none such were available for the present effort. 



III. RESULTS 

A. Total Error in Calculated Mixing Ratio 

The calculated mixing ratio, X , was computed from the usual formula, 

X„ = 0.622 e s(T) 

P - e s (T) 

where e s (T) is the saturation vapor pressure at the temperature T .of the 
saturator, and p is the total gas pressure in the saturator. The values 
of e s (T) were taken from the Smithsonian Meteorological Tables, using the 
temperature indicated by the quartz thermometer, as discussed in Section 
II. B. 2. The total pressure measurement has been discussed in Section 
II.B.3. These two measurement errors were treated as random and independent 
and were propagated according to the usual "root sum-of-squares", giving a 
total percentage error estimate for the calculated mixing ratio of +0,4%. 

B. Total Error in Measured Mixing Ratio 

The measured mixing ratio was simply the mass of water collected, 
divided by the mass of ' dry air expanded. The former, however, had to be 
corrected as discussed in Section II. C. 2. The latter had to be corrected as 
discussed in Section II.A.l. The uncertainty in the measurements of 
collected water were about +0*0? g, from canister weighing procedures, and 
about +0.2 g, from the bag sample check procecure. . The root sum-of-squares 
of these two errors is 0.21 g', or about 1.6% of a total collected water mass 
of 13 g. 

The dry air mass measurement had very small error by comparison. In 
the worst cases, those involving the use of the locally-supplied compressed 
dry air, the total uncertainty in the measurement was about .04 g, from the 
error in the estimate of the original mixing ratio of that air, and about 1 g. 



due to the Ohaus 111 9D balance. The root sum-of-squares of these values 
would be very close to 1 g, giving a 0.11% error estimate for a dry air 
mass measurement of 900 g. 

In the case of a division operation, such as dividing the dry air 
mass into the collected water mass, the percentage errors of the divisor 
and dividend are propagated by a root sum-of-squares rule if they are 
independent and random, giving for the total percentage error of the mea-- 
sured mixing ratio the value 1.6%. As previously discussed, this error esti- 
mate is probably rather conservative due to the conservative way in which its 
major contributor, the bag sample procedure error, was estimated. 

C. Discussion 

Results are given in Table 1. Comments which generally summarize 
important problems or changes are given in Table 2. Note that Runs 1, 4, 
and 5 had problems so serious as to completely eliminate their data from 
consideration. Also, it was suspected in Runs 6 and 7 that the upstream 
canisters were leaking (though it was never possible to detect the leak 
among the many junctions and seals involved); eliminating those canisters 
led to the rather good agreement of Run 8. 

Since the calculated mixing ratio could be given to 0.4%, and the 
measured value to only 1.6%, Runs 2, 3 and 8 would indicate that the satura- 
tor is performing without detectable -errors. Furthermore, the agreement in 
Runs 6 and 7 is good enough, given both the known and suspected errors, that 
no malfunction of the saturator would be indicated. 

The most obvious ways to improve this calibration technique have 
already been discussed; namely, improvement of the canister design and weigh- 
ing procedure (provide, for example, one large canister and a top-loading 
balance of high enough capacity that the trapped water can be weighed with two 
instead of six operations), and improvement of the bag sample check procedure. 



TABLE 1 



Run No. 



Saturator 

Calculated 

Mean Pmb 

Mixing Ratio 



Weight of 
Dry Air, g 



5 x Water Measured Mixing 
n Bag, g Ratio 














TABLE 2 


COMMENTS 


Moisture in bag probably influenced by previous use on room air. No flushing 
flow. This run mainly to familiarize workers with procedures. 

Measured mixing ratio 0.8% higher than predicted. No flushing flow. 

Five minute delay between flushing flow and calibration flow. No detectable 
disagreement between measured, predicted mixing ratio. 

Sump allowed to nearly dry out. 

Added upstream dessi cant .canisters - failure of seals at beginning. 

Bag seal accidentally failed. Calculated mixing ratio 2,8%/higher than 
measured. 

Calculated mixing ratio 2.7% higher than measured, 

Removed upstream dessicant cainisters. Calculated mixing ratio 0.7% higher 


than measured. 




1. Air bottle in place on Ohaus 1119D balance. 



Filled Canister 


Open Saturator 
Covered with 
Insulation 


Quartz Thermometer 
Probe Atop Carbon 
Air Filter Pro- 
vided with Satura- 
tor. 




2. General view of saturator, filter canister, and 
quartz thermometer probe. 




3. Drierite canisters and flowmeter. 
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Drieri te canisters connected to saturator 


NOTES ON SATURATION DESIGN - I: General Considerations 

Dr. P. Squires 

1 . INTRODUCTION 

The purpose of the saturator is to produce a flow of 
order a liter a second with an accurately known mixing 
ratio , by exposing the air sufficiently to wet surfaces at a 
known temperature r 

2. MAJOR DESIGN' ASPECTS 

A. ' Plate Spacing (h) 

If the air is to be exposed for (say) 10 time constants, 
the total volume of the channels is about 5 Fh^, independent 
of the number of such channels (n) . Thus it is desirable to 
keep h small. The volume of the plates is *10 NFht where N is 
the number of air channels and t is the plate thickness. It 
-is therefore desirable also to keep u small (preferably/ N « 1) . 

(see pg. 1 of "Notes on the Saturator Design - 3", in this 
Appendix, "Part 3, for this, derivation) . 

B. Plate Contamination 

Contamination of the wet surfaces by aerosol particles 
does not appear to be a significant factor, but if surfactant 
materials are present and accumulate in- the critical region , they 

ma y cause aexe ter ious. effects, e.g. reducing the rate at which water 
vapor molecules pass through the surface. 

With a design in which n. and h are small (in order to econo- 
mize volume) , the drag of the air flow on the surface of the 
water In the plates is so large that this surface will move down- 
stream with the air, if it is free to do so. This would appear 
to indicate that either the capillary surface of the plates should 
be. in the form of transverse grooves, or else the supply of water 
should be introduced at the upstream end and some surplus water 



drained off at the downstream end, beyond the critical region. 

If transverse grooves are used, and the transverse- water 
flow is expected to transport the surfactants to one side of 
the plates (beyond the effective working humidification region) - , 
the question remains — how rapidly will the surfactant layer 
accumulate- and build back. into the workiricr reaion?' At. present/ 
seem s ke no way to design such a system in such a way 
ensure the required performance , unless the surfactant 
layer is continuously removed on the "far" side of the plates. 

Thus either of the possible designs seems to require that 
so much water be supplied to the plates that there will be an 
overflow to be drained away (either at the "far" side, or at 
the downstream end) , and that this water be drained off in such 
a way that the surfactant layer is carried with it . This latter 
requirement -may prove troublesome: it seems not too likely that 

a capillary wick surface into which the water could drain would 
readily accept the surfactant layer. One method which probably 
would achieve the desired result would be to extract the water 
from the ends of the grooves by means of a. rapidly-moving 
air stream,, since the stress would' be applied directlv-to the 
surfactant layer . This however would' - represent, a very consider- 
able complication of the overall design.. Furthermore, at the 
downstream end, with closely spaced plates, it is doubtful if there 
would be space available to install the required "extraction" 
devices, unless the grooves could be continued around a corner 
on the plate so that the groove ends were of order a centimeter 
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from the saturator air flow. In the case of transverse grooves, 
there probably would be more space available for the extraction 
operation, along the "far" side of the saturator. 

A more radical approach would be to keep the plates wet 
by condensing water- vapor onto them from the air stream. It 
seems reasonable to assume - that' by the use of' activated charcoal 
filters, the long— chain -molecules which form surfactant layers 
could be effectively stripped from the air stream. It may be 
necessary to pass the entire saturator flow (F ^ 10 3 cm 3 sec _1 ) 
through such filters so that surfactant typei molecules originat— 
^- n 9’ ^- n various parts of the system do not accumulate in the 
recycling air. If this is practicable, the only remaining source 
of surfactant would be the water supply to the saturator. If the 
water surface on the plates were being continuously replaced by 
condensation, while the surplus water were being drained away, 
the risk of contamination would seem to be lessened. If the 
contamination remained small enough, it would- be possible to make 
use of capillary wicks to -drain the plates , to a capillary reser- 
voir similar to those which have been planned as sources of water. 
With, this; design, it would seem to be .necessary . to use. a wicking. 
surface consisting., of longitudinal, grooves, since condensation 
cou ^ not be relied upon to keep the plates thoroughly wet in the 
critical downstream region, where the air stream is very close to 
saturation. The capillary drain would then be located at the 
downstream end. 
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It is clear that the concern about the effects of surfactant 
materials is a major driver in the overall design. Provided that 
arrangements can be made to start with the plates thoroughly clean, 
and that evaporation is never allowed to occur from the plates, 
even when the plate temperatures are being, changed, the continued 
renewal of the water, surface in a "condensing saturator" would 
appear to be attractive . Its practicability depends on the design 
of a "vapor generator" which would add water vapor to the air 
stream to raise its dew point above the main plate temperature. 

This generator would need to be flexible enough to be able, to cope 
with transition situations; for example, if the plates were being 
warmed from 5°C to 20°C, unless water vapor could be added to the 
circulating air fairly rapidly, evaporation would occur from the 
plates, causing a reverse flow from the capillary reservoir, and 
increasing the risk of surfactant contamination. 

In connection with a condensing saturator design, it would 
probably be possible to design this in such a way that, with 
-relatively minor changes, it could be operated as an evaporating 
saturator,, using the capillary reservoir as a source instead of 
ns- a .sink. 

C. The Main Plates 

As discussed above-, a single air channel design with closely 
spaced main plates is desirable, since it economizes volume. 

However, as the plate spacing (h) is reduced,, the pressure 
drop in the saturator increases, as does the frictional drag of 
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the air on the water surfaces. This could result in difficulty in 
ensuring the proper irrigation of the* plates. Furthermore, too 
large a- gradient of pressure in the air stream will result in an 
error in the mixing - ratio, which depends on the total gas pressure 
as well as on the partial pressure of water vapor.. 

Some economy of volume may be achieved if the plates, instead 
of being, parallel, .diverge slightly-, downstream. In this way,, it 

may be possible to some extent to combine the advantages of close 
plate spacing (at the upstream, non-critical end) with satisfactorily 
low values of the air velocity and pressure gradient, (in the cri- 
tical region at the downstream end) . If it is desired to mechani- 
cally reinforce the plates by placing . space r s^be tween them, it 
would probably be necessary in this case to make these in the form 
of tapered rails rather than washers. 

D. The Reheater 

The saturator must include a reheater to raise the temperature • 
of the saturated air stream to a level where it .can be allowed to 
flow to other components of the ACPL without risk of condensation. 

The reheater may be in two sections? one at the end of the main 
plates to raise the temperature by a small amount, and a second 
one- 'downstream to raise the.. temperature, further,' to. any- desired 
level . 

Designs which have been discussed so' far feature a reheater 
consisting of parallel plates maintained at a temperature above 
that of the main plates. An insulating section is used to restrict 
the flow of heat back into the main plates. A more attractive 
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^PPtoach may be to heat the air directly by means of electrically 
heated wires stretched across the stream. This would appear to be 
simpler, and would utilize the air stream itself as an insulator, 
possibly leading to a more compact design. 

3 . CONCLUSIONS 

The choices discussed above depend critically on whether it is • 
possible to design an effective and compact vapor generator and direct 
electrical reheater. 
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NOTES ON THE SATURATOR DESIGN - 2 s The Reheater Channel 


1 . INTRODUCTION- 

The simplest reheater design would appear to be one in 
which the channel between parallel main plates (of depth h cm) 
is continued between insulating walls (possibly metallized to- 
avoid electrostatic effects),' and heat is added to the stream' 
electrically, using- taut transverse wires? some distance beyond 
the heater assembly, the channel opens into a transverse duct 
which collects the whole air flow. The material used for the 
heater channel walls is assumed to have a thermal conductivity (k. } 
about 10 times that of air (i.e. 6xl0 _4 cal cm“ 1 deg" 1 sec“ 1 ) . A 
capillary barrier may be included at the downstream end of the 
main plates.. 

2. HEATER DESIGN 

In order to avoid all risk of thermal modification of aero- 
sol particles, it is desirable that the heater wires not be too 
hot - say less than 10 o C- above the air temperature. It is- assumed 
that, at this stage,., it will be sufficient to raise the air stream 
temperature by AT. — 2°C (further heating may be added in a iess 
critical location downstream if desirable) . A 2°C temperature rise 

y 

will require the dissipation of about 2 watts (1/2 cal sec' 1 '} , since 
the mass flux of air is about 1 g sec 1 *. 



For reasons discussed in Section 5 of "Notes on a Saturator 
Design - 3 ", it is desired that the pressure drop through the 
reheater channel (caused by wall drag and wire drag) should not 
be too large; as shown there, this requirement can be met if a. 
design is chosen in which b ^ 30 cm, h = 0.15 cm. The present 
purpose is to discuss, the design of a. reheater channel of the 
same dimensions,' with parallel walls-; it is assumed here that- 
these walls are effectively insulated from environmental influ- 
ences . 

3. HEATER CHANNEL DESIGN 

A. Upstream of- the Heater 

To design this channel, it is necessary to evaluate the longi- 
tudj.rial flow of the heat in the heater channel walls back to the 
main plates. Beyond the heater, the walls will be heated by the 
air, and by the metal duct which is attached to the downstream end 
•of the walls and will no doubt be close to the actual heater air 
temperature (T+AT) . This question is treated here as a one- 
dimensional problem, as if the thickness of the channel wall were 
small compared with its length. 

Considering first the region upsuream or me-, neater., -let 
the- temperature of the main plates be T, and that of the channel • 
walls opposite the heater assembly, T+ST. Heat will be conducted 
back to the main plates, but will also be removed from the reheater 
channel wall by the air flow as it approaches the heater. Let the 
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total effective cross-sectional area of the material in each 

channel wall be A (= b w, where w is the thickness of the channel 

wall) and the thermal conductivity of the wall material, k . 

w 

Clearly, the thinner the wall the better, from a thermal point of 
view. It may be found advantageous to stiffen the wall structure 
by fabricating it with ribs. If so, the cross-section of longi- 
tudinal ribs- should be included in A. _ It is assumed that the 
insulation material (with a thermal conductivity about equal to 
that of air) used to isolate the channel walls from environmen- 
tal influences will not contribute significantly to the longitudi 

nal transfer of heat in the walls. 

Let the origin of x be at the upstream end of the channel 
wall, i.e., at the end of the main plate. The solution due to 
Polhausen for heat transfer from a flat .-plate at temperature T w 
(e.g., Goldstein, vol. 2, 1965, p. 627) gives an expression for 
the total heat transfer: 

1 

Q = 0.686 a 3 k b (T -T) 

cl W 

where' b is the breadth of the plate, T the air temperature at 

infinity, x the plate length, k the thermal conductivity of air 

a 

and u Q the air v^T^^ity of' the free stream. 

Differentiating, one derives an expression- for 'the' rate of 
transfer per unit length at the downstream end: 



dQ = 
dx 


0.343 


1 

a 3 k b 
a 


<V T) 
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This expression may be used to give an- estimate for the 

of heat from the channel walls. Xt tends to overestimate 
hh e local rate of heat transfer because the gas stream is of 
finite depth only . On the other hand , it tends to underestimate 
the local heat transfer because the gas-wall temperature contrast in- 
creases downstream. Further studies are needed, to more clearly define 
the nature of this approximate estimate-. (Underestimating heat 
loss' to the air is conservative, because it will result in 
exaggerating the; heat flux to the main plate. As discussed in 
Section B below, downstream of the heater assembly, the two errors 
cooperate to ensure that the treatment is conservative) , 

Considering conduction in the material of the wall. 


d^T 

dx^ 


0.343 


1 

O 3 k b' (T -T) 
a w 



This may be written 


where y = T^-T and a 



_1 

ay x~ 2 


1 

0.343 a 3 



k b 
a. 



Taking u Q — "the velocity in the' center of the. channel', u = 

* ' -i ^ 

333 cm sec and wxtn k w - lo k a , w .= 1 cm (A = 30 cm 2 ) , it 
results that a - 1.4-. 

The boundary conditions are: y = 0 at x = 0, y = 6T at 

x = where is the distance' from the main plate to the heater. 
Considering the boundary condition at x = 0 , it is found conven- 
ient to write a solution in the form: 
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• 1 2 

3 J 

y * C x a x f(z) = C z f(z) 

_3- 

whcre is an arbitrary constant, 2 = ax 2 , and f(z) is a function 
which may be expanded as 


f (z) = 1 + a, z + 


n n 


+ a 2 z + a z + 


Then, 


2 5 3_ 

y = f(z)+f C,_'a 3 x 2 f(z) 


2 'l- 


y = 


5 x - £ n 

f'<z.) +. | c i “ 3 x 2 f " l2) 


so that the required solution for f (z) is that for: . 

9z f" (z) + 15 f'(z) — 4f (z) = 0 

wh^rp 

f(z) = 1 + a, z + a_ 2 2 ... + a z 11 + .... 
x z n 

and 

f (z) = a 1 + 2a 2 z + 3a 3 z 2 + ... (n+1) a n+1 z n + ... 
f M (z) » (2x1) a 2 + (3x2) a 3 z + n (n+1) a n+1 z n_1 + 


Thus, the coefficients a^ may be found, as follows: 

15 a 1 » 4, 

9(2x1) a 2 + 15x2 a 2 = 4a^ 
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and in general 



Evaluating the coefficients in turn gives: 

a 1 = 0.26667, a 2 = 2.2222xl0“ 2 , a 3 « 8.9787xl0~ 4 
a 4 = 2.1378xl0 -5 , a $ = 3.3534xl0~ 7 , a g = 3.726x10 


The value of C t -. is determined by evaluating the solution at 

3 . 

x = z = z 1 = aA 1 2 (where y = <5T) ? values are given for various 

3 

-T 

values of a ^ in Table X- 



TABLE 

_I 

C t /6T « 

= a *, 3 / 2 

f (z T ) 

*, 2 / 3 f<V> 

1 

2/3 ~ r 

Z , 

1 

1.28981 

1.2898] 

0.77531 

2 

1.62976 

2.5871 

0.38654 

3 

2.0261 

4.2144 

0.23728 

5 

3.0156 

8.8177 

. 0.11341 

10 

'7 ..0378 

'32.667 

0..03061 


The. value of the first derivative of y with respect to 
x at x = 0, is Cj a- 2/3 ,so that the flux of heat (H) from the heatei 
channel wall into the main plate is: 
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2 . Ak a 3 6T 

H = A Vr “ —t 

z l 3 f(z l> 


u o 


Thus, for k w = 6x10 * cal cm ^sec ^deg A = bw — 30 cm 3 , 

|5t? = 333 cm sec’ 1 (a = 1.4,) <$T = 1, this flux is 
2DJtl 


-2 


2.3x10 
2 

Z 1 f ^ z l^ 


For example, if z^ = 3, this amounts to 5.4x10 3 cal sec 1 . 
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B. Downstream of the Heater 

Xn the section of the heater channel downstream of the 
heater wires, a flux of heat passes from the air (now raised to 
temperature T + AT) to the walls. Another source of heat will 
be the transverse duct used to collect the air from the whole 
width of the saturator, which will be attached to the downstream 
ends' of the heater- channel walls.’ This duct is expected to be 
close .to T + AT; and it will be assumed that the outermost ends 
of the channel walls are maintained at this temperature by con- 
duction from the duct. 

As noted earlier, the approximate procedure of differen- 
tiating the Polhausen solution for a plate at constant tempera- 
ture may or may not result in an underestimate of the local, rate of 
heat transfer from the wall to the air in the region upstream of the 
heater, which would ensure that the treatment is conservative. EEow- 
sver, downstream of the heater,- it certainly results in an over- 
estimate of the heat transfer from the air to the wall both because 
the gas is of finite depth and because the gas-wall contrast in tem- 
perature decreases downstream; this is conservative for. the present 
purpose, which is to estimate, the flux of heat to the main elate. 

It is -convenient -in. this instance to locate the origin, of " 
x at the plane of the heater, with x increasing downstream. The 
heat transfer to the walls -from the air is then given by: 
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1 

f§ = 0,343 ff3 kb (T + AT - T )\ 

a . w'V vx 

This influx of heat must be balanced by an increase in the flux 
carried by conduction within the channel wall towards negative x 
(i.e. upstream) . 

Thus : 

d 2 T -I 

— = a (T -T-AT) x 2 
a * w 


and writing 


Y.= T w -T-AT, 


d 2 Y 

dx 2 


= aY x 


1 

’2 


Three boundary conditions designated below as (a) , (b) , (c) 
must be satisfied /.'in the manner described below: 

(a) Y = ST-AT at x =0 

(b) Y ~ 0 at x - £ 2 where is the length, of the channel 

downstream of the' heater. 

... ,dY, ,dy, 

c ' ; l dx' ' “ -'a?)- / to* ensure that the gradient of T 

x=0 ax ,x w 

is continuous at the plane of the heater. (Note that' the origin 
x is different on the two sides of this equation) . 

The solution utilized in the previous section vanishes at the 
origin, and therefore cannot meet condition (a); it must be generalized 
by the addition of an independent solution. The third degree of free- 
dom is then handled by taking account of the fact that physically, 

ST is not an independent boundary "condition , but depends on AT. 
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The required independent solution may be found in the form: 

Y = g (z) 

3- 

where z = ax^ 

Thus-, it is required^ that 

Y = | ax 2 g ( z ) 

JL 

11 3 2 1 9 2 w 

Y- = j ax ..g (z.) + j a x g (z) 

It is therefore necessary to solve: 

» » 

9z g (z) + 3 g (z) - 4 g(z) = 0 

Writing: 

g(z) = 1 + b.z + b_z 2 + .. b z 11 + .... 

±4 n 

g‘ (z) = b ± + 2b 2 z + 3b 3 z 2 +- . . . + (n+1) b n+1 z n + ... 

g (z) = (2x1) b 2 + (3x2)'b 3 z + ... +n (n+1) b ^.z 11 ” 1 + ... f 

it is seen that: 

3b L = - 4 

9-.(2xl);b 2 + 3 ( 2 ) b 2 = 4bj_, etc-. 

and in general, for n > 1, 

9n(n4.1)b n+1 + 3(n+l)b n+1 = 4b n , 

so that b = 4b n 

3 (3n+l) (n+1) 
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For example, b 2 = ^ b l' b 3 = ST b 2 

4 120 b 3 ' b 5 195 b 4, b 6 - 288 b 5 ' b 7 = 399 kg 

Evaluating- the coefficients in turn gives: 

b 1 = 1.3333, b 2 - 0.22222, b 3 = I.4109xl0~ 2 

i > 4 = 4.7031xl0~ 4 , b 5 =■ 9 . 6474xl0 -6 , b fi = I.3399xl0 ~ 7 

b ? = 1.3433xl0 ~ 9 

Thus a general solution is : 

2 

Y = C 2 Z 3 ‘ f (z) + C 3 g(z) 



= C L h.(. 2 ^) 

where h.(z}. = f (z)- +• % zf’.(z) 
Hence , 

h (z, ) <ST 

c = x 

2 2 

2 l 3 f(a 1 ) 
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From condition (b) , 

2 

C, = " z 2 3 f <z 2> C 2 

g Tip 

-z 2 3 f(z 2 ) h(z 1 ) ST 

=• — j 

Z l 3 f(z l ) g(z 2 J - 
From condition (a) 

C 3 = ST -AT' 

Hence , 2 

AT _ t j_ Z 2 f (z 2 } ‘ K(z l ] 

ST ± 2 

Z l 3 f(z l ) 9( z 2 ) 

Consequently , the heat flux to the main plate (H) is: 


H = — £ 


Ak a ST 


z x 3 f (z 1 ) 


a 9 ( z 2 ) AT 


1 r 

z 1 3f(z 1 )g(z 2 )+z 2 3 f(z 2 )h(z 1 ) 


Ak a 
w 


AT 


f 

-z. 


3 • 


L.. f(z l)- t 2 2 f (z 2 ) M z i> 


"gTzTT 


It will be seen that, if z 2 =0, so that ST = AT, this result 
is consistent with that derived at the end paragraph A above. (H = 


2 2 . 

3 3 

Ak^ a ST/ (z 1 f (z 1 ) ) . The second term in the denominator repre- 
sents the beneficial effect in reducing H which is gained by con- 
tinuing the heater channel walls beyond the heater. 
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It is easily shown that 


h (z) = f (z) +; ~ zf ’ (z) 


= 1 + f a i 2 + § a 2 z2 + • • • + — 3n *- -' a n zn + • • • 


where the coefficients of z 11 in h(z) are 


n 


(3n+2) 

2 a n 


1 

2 

3 

4 

5 

6 


6. 6667x1 0'* 1 
8.8889xl0* 2 
4.9383xl0* 3 
1.4965xl0* 4 
2 . 8504x10*® 
3.726xl0* 8 


Values of h(z) 


are given in Table 2 . 


TABLE 2 


z 

. h(z) 

1 

1.7606 

2 

2.7309- 

-3 

3.9462 

-5 

7.2759 

TO:- 

23..-313' 


2 

Values of g(z)- and of z 3 f'(z)/g(z) are given in Table 3. 
(The latter function appears in the expression for H) . 
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TABLE 3 



1 

2 

3 

5 

10 


2.5702 

4.6763 

7.4215 

15.312 

56.480 


0.50183 

0.55324 

0.56786 

0.57586 

0.57838 


From" Table. 3, it can be seen that provided z^ is large 
enough (say, z 2 > 2) , H is essentially independent of .z^ and one 


Ak a AT 

may write H = — — , 

D(z, ) 


where D^) = ^ f (z^ 


z 2 f < z 2> 


Z 2 cor 3 *® s P OI1< is to the situation where the effect of the heat 
flux from the heated air into the wall downstream of the heater 
dominates that due to heat conducted from the transverse duct 
to the downstream encLs of the channel walls - which has been assumed 


here to 
(T- +. AT) 


, # P 

maintain the downstream ends (x = £-) at a temperature of 

^ 2 

) . Substituting a constant value of 0.57 for ^3 f(z 2 )/g (z 2 ) 


2 

yields approximate values of H/(Ak r a 2 3 AT) as a function of z., , 
alone as shown.-, in.. Table 4 ., 
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TABLE 4 


Approximate values of 


H 


Ak^ a AT 


DTTT for Z 2 > 2 


Z 1 


D(z 1 ) 


1 

2 

3 

5 

10 


0.436 

0.241 

0.155 

0.077- 

0.022 


Now, since u = 3F/2bh, 


1 2 


Ak a 3 = (Ak b) 3 (0.343 cr 3 k ) 3 flUA 
w w a \2hvj 

Substituting cr= 0.73, k = 6xl0~ 5 cal cm~ 1 sec“ 1 deg" 1 , 


3 3-1 

F = 10 cm sec , h = 0.15 cm, v= 0.16, 


2 1 
Ak^ a 3 = 2.78xlO~ 2 (Ak„b) 3 


w 


Writing A — bw, where w is the thickness of the heater channel walls 

2 1 


H = 


2.78xl0“ 2 b 3 (w k ) 3 AT 

w 


D(z l ] 


For w = 1 cm, b = 30 cm, k w = 6xl0~ 4 cal cm" 1 sec _1 deg~ 1 this gives: 


H = 


2 . 26xlQ"* 2 AT 
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It is clear from Table 3 that there is little or nothing 
to be gained by increasing above 3. With the design parameter 
values mentioned above (for which a “ 1.4) this implies ~ 2 cm. 
Thus, if an overall channel length of 6 cm is- considered accep- 
table, may be chosen as . large as. 4 cm,, so that with a~ 1.4, 

z 1 % 11, so that % 0.02. For AT = 2°C, this would yield 

-4 - . — 1 1 

H = 9x10 cal sec 

Even is none of this flux were removed by the circulating 

water within the main plate, but the whole flux to the two plates 
-3 -1 

(2x10 cal sec ) were used, to heat the air stream and to evaporate 

additional water into the stream, the .resulting increase in the 

mixing ratio of the air would be only about 2x10 6 g g ^ (about half of 

the heat would be used for evaporation and half for heating the air 

stream) . At the minimum temperature (5 C) , where the saturation 

-3 —1 

mixing ratio at 1 atmosphere is about 5x10 g g , this would result 
in an increase in the mixing ratio of only 0.04%, which is negligible. 

(In fact, the error will probably be much less than this estimate, 
most of the heat. flux from the channel wall being removed by the 
circulating water within the- main plate) . 
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4. 


OTHER PERTURBING EFFECTS 
A. Radiation 


The radiative effects will be overestimated if both the 
heater wires and the main plates are treated as black bodies. 

If they were both at temperature T, the radiative- flux passing 
through the- entrance to the heater channel from the- wires; to the 
'main plates would be of order ( 3 ^-). (mrbd)aT 4 where n is the • 
number of wires and d is their diamtefs; the flux returning from 
the main plate to the wires would be the same. When therefore 
the wires are hotter than the main plates, being at temperature 
T+ST, the net flux from the wires to the plates is of order 
{2crT nhbd ST/l^ cal sec” 1 . 

With the design parameters discussed in Notes on the 

Saturator Design -3, {n * 5, h * 0.15 cm, b = 30 cm, d = 1.3xl0 2 cm, 

ST = 6 °C, ^ = 4 cm ) ' this flux amounts to about 3xl0” 5 cal sec” 1 . 

This is equivalent to less than 0.02 times that estimated as due 

to conduction from the two heater channel walls, and is therefore 
negligible. 

Radiative- heating of the two channel walls by the heater 
wires, treating both as black bodies, is of order 4ntbdcrT 2 6.T cal 
sec which,, -for the values quoted above, amounts to 5 x 1 0 -” 3 cal sec” 1 , 
at the most? a rough estimate can be' made of the resulting error as 
follows. Most of the radiant heat will reach the walls close to the 
heater wires. If it is assumed that all of it is deposited directly 
opposite the heater, the heat flux so introduced may be compared 
with that flowing through the wall towards each main plate of this 
point (H r ), according to the analysis given in Section 3 above. 
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Thus , 


*k - ‘i’x . o = a2/3h( ^ 


whereas = CL a 2//3 

1 


Thus, H — h(z^) H. With the design discussed above, = 
II, and h(z^) ~ 25 (Table 2). Hence (in each wall) . H- ' ~ 2.5xl0~ 2 
cal sec* , whereas- the additional flux deposited by radiation is 
about 2.5x10 cal sec It would therefore appear that radiative 
heating of the channel walls would probably result in an increase • 
m H by about 10%. with the saturator operating at 5C, this- could 
cause the mixing ratio of the air to increase by a further 0.004% 
only. These estimates are likely to be quite exaggerated, since 
they are based on the assumption that none of the heat leaking 
from the channel walls to the main plates will be carried away 
b y the temperauure controlling water flow within the latter) . 

B 

Upstream Diffusion of Heat in the Air 

This may be neglected provided P r > io 2 . with the design 
parameters discussed above/ this nnnrH tion is met: 


R ' P _ ^ 

e r. ' bhv ic 
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s: 


CONCLUSIONS 


The proposed reheater design, with "5 transverse wires 
raised 6C above the temperature of the air stream, has walls 
® cni long, and dispenses with the need for a temperature 
control system for the reheater walls . The total energy dissi- 
pation is about 2 watts ,. and the leakage of hea-t back to the 
niai-n.. plates is trivial; the insulating properties" of the air ' 
flow itself are helpful in this regard. 

As discussed in Notes on the Saturator Design - 3, the 
pressure drop- through the reheater channel (b = 30 cm, h = 0.15 cm) 
is acceptable. 
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APPENDIX II 

3 i MAJOR REQUIREMENTS 



NOTES ON THE SATURATOR DESIGN — 3: Major RGcjuixaments 


1 • Introduction 

Tiis purpose of this section of these notes is to explore major 
features of the design of the saturator, assuming that the longitudi- 
nal gradient of temperature in the saturaror plates is small enough 
to be neglected. . Preliminary investigations having indicated 
that the. advantages gained by using slightly diverging plates 
are not sufficiently great to compensate for the added complexity 
of construction, the discussion assumes parallel plates. The 
Reynolds number of the flow between parallel main plates is 
j— j;, where F (cm sec ) is the flux of air (10 3 cm 3 sec” 1 ) and b is 
the breadth of the main plates. Exploratory calculations indicate 
that the overall length of a single— channel saturator will probably 

order 50 cm; it would appear convenient therefore for packag- 
ing purposes if b were somewhat less. Provided b < 20 cm, R^ < 350, 

and the flow will be laminar. Also, provided b < 50 cm, r p = L 
2 e r bK 

will exceed 10 , so that longitudinal diffusion of heat ‘may be 
neglected. Similarly, if b < 40 cm, the analogous parameter 

exceed 10 , and longitudinal diffusion of water vapor may also 
be neglected. . It will be assumed in- the following discussion that 
a value of b will be' chosen lying between 20 and. 40 cm. It is also 
assumed that the plate spacing (h) will be not less than 0.1 cm. 

Using spacers between the plates to help maintain uniformity of 
spacing, this would seem to approach the practical limit beyond 
which variations in the plate spacing might cause significant chan- 
nelling of the air flow, with some degradation of the saturator 
performance. 



In these circumstances, the approach of the temperature and 

vapor density to equilibrium with the main plates may (far enough 

downstream) be estimated by means of an exponential formula in 

which the residence time is calculated from the mean velocity 

v — F/bh, and the time constant is = h /it k for temperature, 

and x - h 2 /ir 2 D for vapor density. Since D ~ 0.24, k ~ 0.20, 
w 

t > t , and the saturator design may be discussed conservatively 
h w 

in terms of x = x^. 

If the air is to spend n{~10) time constants between the main 

plates, their length L must be given by: 

T _ nhF ~ 

L 5 

ir^tcb 

Thus the area of the plates is proportional to h. From 
this point of view therefore, it is advantageous if h can be kept 

small . 


2. The Direct Effects of Pressure Gradient on the Mixing Ratio 
The longitudinal gradient of pressure increases rapidly as 

h decreases: 

d£ = -12r|F 
dx . bh 3 

This pressure fall downstream affects the value of the mixing 
ratio (w) . If the vapor pressure is in equilibrium with the channel 
walls, as a result of the downstream pressure fall alone. 


dw 

dx 


12r|Fw 

bh 3 p~ 


2 



At the downstream end, therefore, the mixing ratio would be 
elevated above its equilibrium value by 6w, where 


Sw _ xv dw 
w w dx 


!2nF 2 

TT 2 icb 2 h 2 p 


Substituting r\ = 1.6xl0~ 4 (cgs) , P = 10 3 cm 3 sec _1 , k = 0.20, 
p .= 10 6 dynes . cm" 2 , 


— = 9*7x10 
W ~h Z h 2 

Provided b > 20 cm, h > 0.1 cm, the proportional error in w 
caused directly by the downstream fall of pressure will be less 
than 0.03%,. which is negligible. 


3-. The Irrigation of the Main Plates 
A. Using Longitudinal Grooves 


If the wicking surface consists of square longitudinal grooves 
of side s, spaced with their centers 2 s apart, the longitudinal 
forces acting on the strip of water in- a gr,oove are those due to 
- (,a) ; • the .downstream . decrease of air pressure, i.e. 


-Ls 2 dp ; I2nF 2 rf s 2 

dx 2 ,22 
tt Kb h 

(b) the drag of the air on the water surface Lsri (^) 


6FsnL = 
w bh 2 


6nsOF 

~2 — -j- , and (c) the drag exerted by the walls of the groove on the 

it Kb h 


3 


water moving' through it. Force (c) is small compared with (b) ; the 

total, flow of water in the grooves , (F ) is expected to be of order 

10 g sec at the most. . This flow is divided between grooves. 

The velocity gradient at the groove wall is therefore of order 
12F W 

— .- The total force acting on the water in a groove due to wall 
bs 

• 36 F n L 

drag on three sides is- therefore of order — , w — . The ratio of- 

. bs 

force (c) to force (b) is* therefore 


* F if' t.2 
w w h 

f n T * 
s 

Provided h j 0.3 cm, s { 2.5xl0~ 2 cm, this ratio is less than 10% v 

and force (c) may be neglected. (in fact, as will be seen later, 
there seems no reason why h should exceed 0.15 cm). 

Forces (a) and (b) (the affect of which are additive) must 

be dominated by the capillary* stress exerted ;by the groove*, which 

4cr 

can exert a pressure of — corresponding to a force on the strip 
'of water of 4as. 


Thus, the condition to be met is: 


6nF 2 r) , 2s s _ 

. 2 '2. (1+ E~* ) 4(7 

tt icb- h 


Thus 


h * 

, provided st^-, as seems certain to be the case, the actual 


magnitude of s is not of dominant importance in this connection, 
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If then the maximum capillary stress which can be exerted by the 
grooves is to dominate the disturbing influence by a factor of 5 , 
it is required that 


b 2 h > -f— n • 170 

^ 'V' 

ir i<0 

for n = 10, F = 10 2 cm 2 sec~^. 

' This requirement specifies minimum values of h for any value 
of b, and consequently, determines also minimum values of L- 


Tt Kb 


. Values of h . 

mm 

and L . are given as 
mxn 3 

a function 

(n = 10, F =• 10^ cm 2 sec~^- for 20 < b < 

4Q . 


TABLE 1 


Possible Designs Using Grooves with 

. n 

s< x 

■ b 

h . 
mm 

L min 

(cm) 

(cm) 

(cm) 

20 

0.43 

109 

25 

0.27 

55 

30 

0.19 

32 

35 

0.14 

20 

40 

0.11 

14 

Using a Sintered 

Metal Wicking Surface 



'If the wicking surface is a sintered capillary material, 
the drag - exerted by the., air -will be- resisted locally by the sur- 
face. The capillary forces need therefore resist only the pres- 
sure gradient, and the condition to be fulfilled becomes: 


5 



<<4 a 



where e is the -pore diameter of the wicking material. If capillary 
forces are to dominate by a factor of 5, 


, 2 . 2 , „ 2 
k h y k5nF r) ^ 170 
s 2 

- TT ICO" 

-3 

Since e can be as small as 10 cm, this condition places no signifi- 
cant limitations on the design parameters b and h. The restrictions 
on the design indicated in Table I may therefore be avoided by using 
a sintered capillary material instead of grooves; the choice of 
h is then effectively independent of that of b, and the length 

s * nbP 

is given by L = — . With n = 10, this results in the designs 

tt Kb 

shown in Table 2 . 

TABLE 2 

Possible Designs Using a Sintered Metal Wicking Surface 


L (cm) 


b (cm) 

h = 0.1 cm 

h = 0.15 cm 

20 

25 

38 

30 

17 

25. 

40 

13 

19 
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4. Energy Consumption and Dissipation 

The energy dissipated per second in forcing the air flow 
through the main plates is 


T2nF 3 n _ 9. 7x10 6 

FAp “ 2 ,22 ^rr 

IT Kb h- bn 


■ 2 2 

The minimum value of b h occurring in the sample designs 

4 

shown in Tables 1 and 2 is 4 cm . Therefore the maximum energy 

dissipation is 2.4x10^ erg sec” 1 , or 0.24 watt. This heat will 

eventually be deposited in the water flow used to control the 

2 -1 

plate temperatures, which will be of order 10 g sec . The 
resulting temperature rise of the water as it leaves the plates 
will therefore be of order 6xlO”^°C, which is negligible. 

The heating of the air stream by the viscous dissipation of 
energy is relevant only to the degree that it results in an increase 
in the temperature of the water surface, since this could 
perturb the vapor pressure. After' several time constants, the 
energy flux density to each wall from the air amounts to: 

FAp 
SSL; 


wV 9^00* cnf^eo- 1 


2.. 3x10 


, -2 -1 
cal cm sec 
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2-3 5 

The minimum value of b h .in Tables I and II is 0.4 cm 3 ; 

hence the maximum flux density occurring in any of the sample 
. . —5 -2 -1 

desxgns is 6x10 cal cm sec . Even if the thermal resis- 
tance between the wall surface and the temperature controlling 
water flow is equivalent to that of 1 cm. of material of thermal 
conductivity 0.05 cal cm ^ sec ^ , (such as stainless steel), the 
resulting temperature rise of the wall surface is onlylO~ 3 C, 
which is negligible. ' 

Thus, the vi-scous dissipation of energy within the air 
stream does not seem likely to restrict the design. 

5. The Measurement of Pressure 

At the end of the main plates, the vapor pressure (e) should 

be accurately equal to the equilibrium vapor pressure at saturation 
at the plate temperature. In order to deduce the mixing ratio, both 
e and p, the total pressure, must be known at the same point. 

The mean velocity of the air between the main plates is ^ , 
which for the designs shown in' Tables 1 and 2 has a maximum 
value of 500 cm sec , corresponding to a dynamic pressure head 
of only about 0.1 mb. Thus the measurement of pressure at this 
point .does not present any problem .from this point of view. 
Nevertheless-, it would seem desirable .to make the measurement 
downstream of the reheater where the air will be. unsaturated, in 
order to avoid any possibility that the measurement might be ' 
perturbed by water creeping into, or condensing in, the pressure 
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sampling tube. If this is to be done, it is essential that the 
pressure change through the reheater be small, so that the total 
pressure at the end of the main plates can be accurately esti- 
mated . 

Preliminary studies of a rather compact and simple "direct" 
reheater are described i‘n "Notes on the Saturator Design — 2,” It 
consists of electrically heated transverse wires (of length b) , and 
it appears that a total reheater channel length of 10 cm would be 
more sufficient, in such a reheater, pressure drop would be 
caused by drag exerted by the walls, and by the transverse wires. 
The pressure gradient due to the drag exerted by reheater channel 
walls spaced h cm apart is: 

- dp _ 12gF 
dx bh 3 

Since the minimum value of bh 3 in Tables 1 and 2 (b = 20 cm, 

^ — 0.1 cm) is 0.02, the pressure gradient in these cases does 
not exceed 600 r\F . Thus, for a channel 10 cm long, the pressure 
change will not exceed 10 dynes cm~ (1 mb) . Since this 
pressure change is relatively constant, it could probably be 
estimated to an accuracy of order 10%, reducing the error to 
0.1 mb. However', a design such as that in the second line of 
Table 2 , in which b = 30 cm, h = 0.15 cm so that bh^ - 0.10 
and the total pressure fall is only 0.2' mb would appear to be 
preferable. 
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Assuming that the transverse heater wires are located centrally 
and spaced so that they interact independently with the air stream, 
both the heat transfer to the air and the drag depend in known ways 
on the Reynolds number of the flow about the wires (R g ) , as shown 
in Fig. 1 by the uppermost full curves. 

If the acceptable temperature rise of N parallel transverse 
wires is AT, and it is desired to raise the temperature of the 

* *** A » — < XX 

air; stream- by S-T, the Nusseit number (~, where H is the heat 
transfer per unit length per degree per second) must be such 


that; 


H _ 

ST NbJcAT 


Also, the pressure decrease caused by the drag of the wire is; 


6p * 


9NdpF 2 C D 

2™3 
8b z h 


where d is the diameter of the transverse wires. 

Assuming for the moment that these relations can be satis 
fled with a suitably small integral value of N (the number of 
transverse wires) , N may be eliminated to give: 


Sp 


9p 2 ’c p F 3 aST- kC D 


8kh 3 b 3 AT 


H 


In Fig.' 1,. values o’f‘ log. Q ' are pitted, as isolated 


*10 


H 


heavy points, and the (dashed) line; 


log 


10 


C D k * 0.62 
H 


0.88 log 1Q R e 


has been fitted by eye to them in the region 0 < log 10 - e 

C ]c 

with an accuracy of hbout 15% in D 

H 

accuracy , we may write ' 


R < 2, 
e 

Thus , to a sufficient 
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^ and p- for a Cylinder 

(H = heat transfer coefficient per unit length, 
k = thermal conductivity of fluid) 

Sources: Goldstein and Roberts, respectively 



Figure 1 



C D k = 4.2 (R*) °* 88 
— e 


over this Reynolds number interval. 


Thus, 6p = 


4 7 P 2 o p P 3 d6T fbh ^ 

= '* — T^3 Fd 

kh b AT J 


Q.88 


= ‘ 4.7 


p 2 c d°* 12 6Tv° 

P 

kAT 


f— 1 

[hbj 


2.12 


Clearly, ‘to minimize «p, it is somewhat advantageous it a is 

chosen as small as possible; a convenient minimum size is five 

—2 

thousandths of an inch (1.3x10 cm) . 

Then, assuming reasonable values such as AT = 6 C, ST = 2 C, 

2.12 


6p =1.7.5xl0- 4 .[Is] 


In Tables 1 and 2, the minimum value of bh is 2, and hence 

does not' exceed 500 cm sec” 1 .' Hence, Sp < Q.6 mb. In a design 
bh 

with h = 0.15 cm, b = 30 cm, which, as pointed out above, seems 

preferable, <5p ~ 0.1 mb; assuming that this pressure change can be 

estimated to 10%, the residual error (~0;01 mb) is negligible. 

' * —1 

In the case .b = 30 cm, h = 0.15 = ^33 cm sec 

R * = 333x1. 3xl0” 2 _ 2? 


V 


The corresponding- value of ^ is about 10. -Thus, the required 
value of N is: 

Fp c ST k 

w = 2 (£.) 

N kbAT, V 


% 5 


In summary, for h > 0.1 cm, the pressure change caused by 
wire drag does not seem to be a serious problem. However, that 
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caused by wall and wire drag can be somewhat larger than is 

desirable, and it appears advisable to adopt a design such as 

✓ 

b = 30 cm, h = 0.15 cm, leading to a pressure drop of less than 
0.3 mb, which could probably be estimated to within 0.03 mb. 


6 .. Aerosol Behavior 

If the aerosol is injected upstream of the saturator, haze 
droplets will form on the particles as they pass between the main 
plates, so that the total water content of the gas passing out of the 
saturator will slightly exceed the saturation value. 

It is assumed here that the aerosol consists of a soluble salt 
such as NaCl, and that cumulative distribution of critical super- 
saturations is given by N = cS£, the range of S extending from a mini- 

mum value (s) up to values exceeding 10~ 2 (1%). The value of s 
is taken as that corresponding to a maximum dry particle radius 

-5 

r m — 10 cm. For an aerosol of NaCl particles, this implies: 



where A = 


• 32 M M Q 2 ff 3 

27 1 p p Po 2r3t3 


where r is cm, s is in absolute units , and M is the molecular weight 
of NaCl (57) , that of water, '<j the surface tension, p ^ the density 
of the dry particles (2.16 g cm" 3 ), p . that of water, and i the van't 
Hoff factor (taken to be 2) . The value of A depends somewhat on 
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temperature? at 20C, A t 1.4xlO' 22 , and at -20C, A * 2.6xl0” 22 . 
Hence, over- this temperature range, s varies from 3. 7x10 ” 4 (0.04%) 
to 5.1xl0“ 4 (0.05%). . 

If sufficient time were available, the haze droplets formed 
on the aerosol particles would grow close, to their equilibrium 
radii (r Q )-, given bv: 


4 M a 

where B = 2 

3/1 R p w T 


The value of B is also somewhat dependent on temperature; at 20C, 

B 2 -4.2x10 8 , and at -20C, B « 5..2xl0~ 8 . 

The total density of water substance (g cm” 3 of air) which is 
present in the form of haze droplets (p^) would then be: 


P h = | irckp B 3 


m 


s^is 

c c 


s 


■ 4'iTcp. b; 

. O 


(3-k)' 


, (k-3) 


3 k— 3 

since << s for values of' k which are likely to occur (k < l)o 

The term B s ^ is temperature, dependent, varying as 

■ CT |(*-d 

<fr 
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Since -is a decreasing function of T, for- k < 1, p h is 
larger at higher temperatures. It is sufficient therefore to 
evaluate at 20C, where: 


= 3.1x10 


-22 ck 
( 3-k) 


k-3 


For-ari .aerosol, .in which N . = 1 0 \m' 3 at S' = 1-6" 


_ _ « 3+2k , 

c - 10 ,. and 


P h - 3.1X10- 13 jgzy (10 2 s ) k ’ 3 

Substituting s - 3.7x10 4 , p^ is found to be about 10~"' I ' 0 q cm” 3 

if k = 1, and 2.5xl0” 10 g cm -3 if k = 0.5 (at 20C) . 

The density of saturation water vapor at -20C is about 
—7 —3 

9x10 g cm . Thus, even at such a low temperature, the water in 
the form of haze droplets cannot represent an increment in the total 
mixing ratio of as much as 0.03%, which is negligible. Thus, it 
is not necessary to investigate whether the haze droplets have 
time to grow to sizes close to equilibrium sizps *f- 

saturation 

This would not. be the case if the aerosol size distribution 
extended' to much larger sizes than r = 10~ 5 cm. Thus, if the spec- 

tr 3/2^ Xten< ^ e< ^ t0 r P = ^ ^ cm ' s would be decreased by a factor of 
10 , and p^ by a factor of 10 3 for k = 1. In that case, it would 

be important to carefully evaluate the growth of haze droplets 
during their passage through the saturator, and the need to know 
the value of p^ with appropriate accuracy might well place serious 
restrictions on the .saturator design. 
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7 . Entrainment of Liquid Water- by the Flowing Air 

Water will not be entrained from the wet walls if the forces 
exerted by the air flow on an "elementary droplet" or nascent 
wave are small compared with those due to surface tension. Let 
the diameter of the elementary droplet be d^cm. Then (the gradient of 

velocity at the wall, being 6. v/h) , the aerodynamic stress exerted on i 

is of order - ‘2 -2 

6r|vdi .. 6riFdi 

h , ,2 

bh 

The capillary forces are of order cd^ Hence , entrainment 
would not be expected to occur if 


6flFd 

2 

bb7o 


■!<< 1 , 


The minimum value of bh in Tables 1 and 2 is 0.2, so that 

7xl0" 2 d-. 

bh ' a 1 

Since d^is comparable with the dimensions of the characteristic 
elements making up the wall surface, it can be of order 10~ 2 cm, 
so that 

6riFdi . 

' —n~< 10 ^ 

bh, a 

There, seems therefore- no danger that, water droplets could be 
entrained from the wall surface. 
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8 . 


Conclusions 


The overall saturator size depends chiefly on the plate 
spacing, h. For values of h > 0.1 cm, there appears, to be no 
problem arising from the direct effect, of the pressure gradient 
2^' on water vapor mixing ratio; or in connection with the 
irrigation of the plates, provided, a suitable capillary wicking 
material is used; or from the dissipation of energy within the 
air stream. 'However, If longitudinal grooves are used as wick- 
there will be certain limitations on the design in connec- 
tion with satisfactory irrigation of the plates . 

In these circumstances, it seems questionable whether the 
use of longitudinal grooves is advisable; in addition, any sur- 
factant accumulating on the water surface would probably accumu- 
late in the critical downstream portions of the grooves, which 
may be undesirable. If a suitable wicking material is used 
instead, surfactants will probably be immobilized. With such a 
surface, the choice of h(>0.10 cm) seems to be unrestricted, and 
then the overall dimensions may be determined so as to provide 
a reasonable aspect ratio for packaging such as 3:2. If the 

total length of the. entry and reheater sections is L , it would 

a 

then be desired thai 

= I b 


■ Since L = 


nhF' 

2 . 
ir Kb 


, this leads to a quadratic equation, the rele- 


vant solution of which is: 


b 


~ L + 
3 a 



2nhF 

T~2~ 

3ir k 


1/2 
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Another, significant consideration, for the choice of a 
design is that it would appear to be desirable to make the 
measurement of total pressure downstream of the reheater channel. 
As discussed at the end of Section 5, this suggests the choice of 
parameter values b = 30 cm, h -0.15 cm, which correspond to main 
plates of length L. = 25 cm (Table 2) . Assuming, as seems likely, 
that 20 cm, the resulting' aspect ratio of the whole apparatus 

is then 45 : 30 , or 3:2. 

This design therefore appears to meet all the requirements 
discussed. 
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